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Abstract 

We critically examine the QCD predictions for the Q 2 dependence of the electron- 
proton deep-inelastic structure function F2(x,Q 2 ) in the small x region, which is being 
probed at HERA. The standard results based on next-to- leading order Altarelli-Parisi 
evolution are compared with those that follow from the BFKL equation, which corresponds 
to the resummation of the leading log(l/x) terms. The effects of parton screening are also 
quantified. The theoretical predictions are confronted with each other, and with existing 
data from HERA. 
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The first measurements of the proton structure function F 2 (x, Q 2 ) at small x have been made 
by the Hl|l| and ZEUS0 collaborations at HERA. A striking increase of F 2 with decreasing x is 
observed which is consistent with the expectations of perturbative QCD at small x as embodied 
in the BFKL equation ||. This equation effectively performs a leading a s log(l/a;) resummation 
of soft gluon emissions, which results in a small x behaviour F 2 ~ x~ x with A ~ 0.5. 

The data at Q 2 = 15 and 30 GeV 2 are shown in Fig.l, together with a representative set 
of predictions and extrapolations, whose distinguishing features we elucidate below. These 
curves fall into two general categories. The first, category (A), is phenomenological and is 
based on parametric forms extrapolated to small x with Q 2 behaviour governed by the next-to- 
leading order Altarelli-Parisi equations. The parameters are determined by global fits to data 
at larger x (examples are the curves in Fig. 1 labelled MRS(D'_)|4j], MRS(H)|| and, to some 
extent, also GRV||, but see below). The second approach, denoted (B), is, in principle, more 
fundamental. Here perturbative QCD is used in the form of the BFKL equation to evolve to 
small x from known behaviour at larger x (e.g. AKMS0). In other words in approach (A) 
the small x behaviour is input in the parametric forms used for the parton distributions at 
some scale Q 2 = Qq, whereas in (B) an x~ x behaviour at small x is generated dynamically 
with a determined value of A. Of course in the phenomenological approach, (A), it is possible 
to input a BFKL- motivated small x behaviour into the starting distributions (e.g. MRS(D'_) 
and MRS(H) have xg, xq sea ~ x~ x with A = 0.5 and 0.3 respectively). Since the x~ x behaviour, 
for these values of A, is stable to evolution in Q 2 we may anticipate that it will be difficult to 
distinguish approaches (A) and (B). However the Q 2 behaviour (or scaling violations) of F 2 is, 
in principle, different in the two approaches. 

The Altarelli-Parisi Q 2 evolution is controlled by the anomalous dimensions of the splitting 
functions (and by the coefficient functions) which have been computed perturbatively up to 
next-to-leading order. On the other hand the BFKL approach, at small x, corresponds to an 
infinite order resummation of these quantities, keeping only leading log(l/x) terms. Summing 
the leading log(l/x) terms, besides generating an x~ x behaviour, gives its own characteristic Q 2 
dependence. One of our main purposes is to study whether or not the BFKL behaviour, which 
is more theoretically valid at small x, can be distinguished from the approximate Altarelli-Parisi 
parametric forms which neglect the log(l/x) resummation. 

If we were to assume that Altarelli-Parisi evolution is valid at small x then 

8F>M<) 2Eet ^m f^pJAy^Q*) + ... , (!) 
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and hence the Q 2 behaviour of F 2 can be varied by simply exploiting the freedom in the gluon 
distribution at small x. However the situation is much more constrained when the BFKL 
equation is used to determine the (unintegrated) gluon distribution f(x, k 2 ^). Then F 2 may be 
calculated @ using the /^-factorization theorem M 



T~ 

T 



F 2 (x,Q 2 ) = I ^ l^ff(^,k 2 T )Fi \x',k 2 T ,Q 2 ) (2) 

rv'j 



where x/x' and k^ are the longitudinal momentum fraction and transverse momentum that are 
carried by the gluon which dissociates into the qq pair, see Fig. 2. F 2 ^ is the quark box (and 
crossed box) amplitude for gluon- virtual photon fusionf?]]. 
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In order to gain insight into the different possible Q 2 dependences of F 2 it is useful to 
introduce the moment function of the (unintegrated) gluon distribution 



f(n,k 2 T ) = / dxx n - 2 f(x,k 2 T ). 
Jo 

The evolution of the moment function is given by the renormalization group equation 

f(n,k%.) = f(n,kl) exp 



k r dk' 2 

—j(n,a s (k%)) 



k 2 k 



(3) 



(4) 



where the anomalous dimension j(n,a 8 ) is known. From eq.(3) we see that the behaviour 
at small x is controlled by the leading singularity of f{n,kj) in the n plane. In the leading 
log(l/x) approximation 7(71, a s ) is just a function of the single variable a s {k^)/{n — 1) and is 
determined by the BFKL kernel. Its value is such that[B[ 



Sa s (k 2 ) 
7r(n — 1) 







is satisfied, with 



(5) 



(6) 



tf( 7 ) = 2*(l)-*( 7 )-*(l- 7 ), 

where \I/ is the logarithmic derivative of the Euler gamma function. 

For fixed a s the leading singularity of f(n, k 2 -) is a square root branch point at n — 1 + \l 
where \ L = 3a s K(^)/7r = 12a s log2/7r. Comparing with eq.(5) we find that 7(1 + A_l,« s ) = ~. 
Thus, from eq.(4), it directly follows that 



f{x,k 2 T ) ~ {k 2 T )*x 



(7) 



Since jk\ in eq.(2) is simply a function of k^/Q 2 , this leading behaviour feeds through into 
F 2 to give 



F 2 (x,Q 2 )~(Q< 



\2X 



(8) 



where in (7) and (8) we have omitted slowly varying logarithmic factors. 

Formula (4) is valid for running a s , provided n remains to the right of the branch point 
throughout the region of integration, that is provided n > 1 + 12a s (fco)log2/7r. (For smaller 
values of n the k\ dependence of f(n,k^) is more involved [TO]-) For running a s the small x 
behaviour of f(x, k^) is controlled by the leading pole singularity of f(n, k^) which occurs at 
n = 1 + A, where now A has to be calculated numerically |nj. A value of A ~ 0.5 is found, with 
rather little sensitivity to the treatment of the infrared region of the BFKL equation []12]. The 



k\ dependence of / (and hence the Q 2 dependence of F 2 ) is determined by the residue (3 of this 
pole. Using eq.(4) we have 



where 



(3{k 



exp 



4 dk'n 



L/2 



■ 7 (1 + A,a s (^)) 



(9) 



(10) 
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From the above discussion we see that this form is valid provided k\ > > k 2 (A), where k 2 (A) 
satisfies the implicit equation A = 12a s (/t 2 (A))log2/7r. Similarly, provided that Q 2 > k 2 (A), we 
have 

F 2 (x,Q 2 ) ~ [3{Q 2 )x-~\ (11) 

up to slight modifications which result from known Q 2 effects embedded in . Note that for 
Q 2 ~ /t 2 (A) we should again get an approximate (Q 2 ) 5 behaviour of F 2 (x, Q 2 ), although it may 
(at moderately small values of x) be modified by the non-leading contributions. Here we are 
also interested in Q 2 < k 2 (X) and then the form of f3 is more involved |10[ . 



In the leading log(l/x) approximation the anomalous dimension, 7(71, a 8 ), is a power series 
in a s /(n — 1). For the BFKL approach 7(71, a s ) contains the sum of all these terms. If only the 
first term were retained then the Q 2 behaviour would correspond to Altarelli-Parisi evolution 
from a singular x~ x gluon starting distribution with only g — > gg transitions included and with 
the splitting function P gg (z) approximated by its singular 1/z term. 

It is useful to compare the Q 2 dependence of F 2 which results from the theoretically mo- 
tivated BFKL approach, (B), with that of the Altarelli-Parisi Q 2 evolution of approach (A). 
For Altarelli-Parisi evolution the Q 2 behaviour of F 2 depends on the small x behaviour of the 
parton starting distributions. If we assume that the starting distributions are non-singular at 
small x (i.e. xg(x,Ql) and xq sea (x,Ql) approach a constant limit for x —* 0), then the leading 
term, which drives both the Q 2 and x dependence at small x, is of the double logarithmic form 



F 2 (x, Q 2 ) ~ exp 

where 



2{£(Q 2 ,Q 2 )log(l/a;)}: 



Q 2 dq 2 3a s {q 



(12) 



t™><n=U 1 * ■ (13) 

From (12) we see that, as x decreases, F 2 increases faster than any power of log(l/x) but slower 
than any power of x. 

If, on the other hand, the starting gluon and sea quark distributions are assumed to have 
singular behaviour in the small x limit i.e. 

xg(x, Ql),xq se3u (x, Q 2 ) ~ x~ x (14) 

with A > 0, then the structure function F 2 (x, Q 2 ) behaves as 

F 2 (x,Q 2 ) ~ x - x h(Q 2 ) (15) 

where the function h(Q 2 ) is determined by the corresponding anomalous dimensions of the 
moments of the (singlet) parton distributions at n = 1 + A, as well as by the coefficient functions. 

We emphasize again that, in contrast to the BFKL approach, for (next-to-leading order) 
Altarelli-Parisi evolution the relevant quantities which determine h(Q 2 ) are computed from the 
first (two) terms in the perturbative expansion in a s . Thus terms are neglected, which may in 
principle be important at small x, corresponding to the infinite sum of powers of a s / (n — 1) in 
7 (and in the coefficient function). 
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Note that in both cases (i.e. eqs.(12) and (15)) Altarelli-Parisi evolution gives a slope of the 
structure function, dF 2 (x, Q 2 ) / d\og(Q 2 ) , which increases with decreasing x. The MRS(D'_)||] 
and MRS(H)|J extrapolations are examples of (15), with A = 0.5 and 0.3 respectively. On 
the other hand, the behaviour of F 2 obtained from the GRV|| partons is an example of (12). 
In the GRV model the partons are generated from a valence-like input at a very low scale, 
Qq = 0.3GeV 2 (and then the valence is matched to MRS at much higher Q 2 ). Due to the long 
evolution length, £(Qq, Q 2 ), in reaching the Q 2 values corresponding to the small x HERA data 
the GRV prediction tends to the double logarithmic form of (12). The GRV model is probably 
best regarded as a phenomenological way of obtaining steep distributions at a conventional 
input scale, say 4GeV 2 , since the steepness is mainly generated in the very low Q 2 region 



where perturbative QCD is unreliable 11311 . Note, however, that the steepness is specified by the 
evolution and is not a free parameter. In fact, in the region of the HERA data, the GRV form 
mimics an x~ x behaviour with A ~ 0.4, although for smaller x it is less steep. 

To summarize, we have discussed four different ways of generating a steep x behaviour of 
F 2 (x,Q 2 ) at small x, each with its own characteristic Q 2 dependence: the BFKL fixed and 
running a s forms, (8) and (11), the Altarelli-Parisi double leading logarithmic form with a 
long Q 2 evolution, (12), and finally Altarelli-Parisi evolution from a steep x~ x input, (15). 
Examples of such forms are, respectively, the fixed and running a s AKMS predictions^ [P2|| , 
and the GRV|| and MRS(H)|| extrapolations. Their Q 2 dependences are compared with each 
other in Fig. 3 at given values of small x in the HERA regime. For reference the MRS(D'_)|4]] 
extrapolation is also shown. The theoretical curves are calculated either from eq. (2) (where / is 
the complete numerical solution of the BFKL equation obtained as described in ref.|12|) or from 
the full next-to-leading order Altarelli-Parisi evolution. We also show, in Fig. 3, Hl|l4| and 
ZEUS|| measurements of F 2 made during the 1992 HERA run, corresponding to an integrated 
luminosity of 25nb~ x . Only the statistical errors of the data are shown. Measurements will be 
made with much higher luminosity, and at smaller x values, in the future. 

Several features of this plot are noteworthy. First, if we compare the data with the u x~ x 
dependences" of the Altarelli-Parisi forms of MRS(D'_), GRV and MRS(H) (which have re- 
spectively A =0.5, "^0.4", and 0.3), then we see that MRS(D'_) and GRV are disfavoured. So 
we are left with MRS(H), which, in fact, was devised simply to reproduce^] the HERA data of 
refs.il. 

Second, we see that the AKMS prediction (which pre-dated the HERA data) is, like 
MRS(H), in good agreement with the x and Q 2 dependence of the data. In principle, it is 
an absolute perturbative QCD prediction of F 2 (x,Q 2 ) at small x in terms of the known be- 
haviour at larger x, but, in practice, the overall normalization depends on the treatment of 
the infrared region of the BFKL equation]?], |j~2 |. We can therefore normalise the BFKL-based 
predictions so as to approximately describe the data at x — 0.0027 by adjusting a parameter 
which is introduced |12|] in the description of the infrared region. For the running a s AKMS 
calculation, this is achieved if the infrared parameter k 2 w 2GeV 2 (with k 2 = lGeV 2 ), in the 
notation of ref. [12] . Strictly speaking, within the genuine leading log(l/a;) approximation the 
coupling a s should be kept fixedQ. We therefore also solved the BFKL equation with fixed a s , 

1 See also the partons of the CTEQ collaboration which have A— 0.27[]T5|| . 

2 The use of running a s has the advantage that then the BFKL equation reduces to the Altarelli-Parisi 
equation in the double leading logarithm approximation when the transverse momenta of the gluons become 
strongly ordered. 
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choosing a value a s =0.25 so as to have a satisfactory normalization. The resulting Q 2 depen- 
dence of F 2 (x, Q 2 ) turned out to be almost identical to that calculated from the solution of 
the BFKL equation with running a s . For clarity, we therefore have omitted the fixed a s curve 
from Fig. 3. Also a background (or non-BFKL) contribution to F 2 has to be included in the 
AKMS calculation^; this explains why MRS(H), with A=0.3, and AKMS, with A « 0.5, both 
give equally good descriptions of the HERA data. However, by the smallest x value shown we 
see that the BFKL-based AKMS predictions for F 2 begin to lie significantly above those for 
MRS(H), due to this difference in A. 

A third feature of Fig. 3 is the stronger Q 2 dependence of the AKMS predictions as com- 
pared with the MRS and GRV extrapolations which are based on Altarelli-Parisi evolution. 
This we had anticipated, with a growth approaching (Q 2 ) 5 for BFKL as compared with the 
approximately linear logQ 2 behaviour characteristic of Altarelli-Parisi evolution. In reality, at 
the smallest x value shown we find that the AKMS growth is reduced to about (Q 2 ) 3 , due to 
the fact that F 2 (background) is still significant. Although we see that the BFKL and Altarelli- 
Parisi Q 2 behaviours are quite distinctive, to actually distinguish between them will clearly be 
an experimental challenge, particularly since Q 2 ~ 15GeV 2 is the kinematic reach of HERA 
at the lowest x value shown. Recall that the BFKL and Altarelli-Parisi equations effectively 
resum the leading log(l/a;) and log(Q 2 ) contributions respectively. Thus the BFKL equation 
is appropriate in the small x region where a s log(l/x) ~ 1 yet a s \og(Q 2 / Ql) <C 1, where Ql 
is some (sufficiently large) reference scale. If the latter were also ~ 1 then both log(l/x) and 
\og(Q 2 /Ql) have to be treated on an equal footing|T7|], as is done, for instance, in the unified 
equation proposed by Marchesini et al.[[L8[]. For this reason we restrict our study of small x via 
the BFKL equation to the region 5 ~ Q 2 ~ 50GeV 2 . As it happens, the very small x HERA 
data lie well within this limited Q 2 interval. 

So far we have neglected the effects of parton shadowing. If, as is conventionally expected, 
the gluons are spread reasonably uniformly across the proton then we anticipate that the 
effects will be small in the HERA regime fl2||. For illustration we have therefore shown the 



effects of (speculative) "hot-spot" shadowing, corresponding to concentrations of gluons in 
small hot-spots of transverse area ■nR 2 inside the proton with, say, R = 2GeV _1 . In this case, 
to normalise the predictions at x=0.0027, we need to take the infrared parameter k 2 ~ 1.5GeV 2 . 
With decreasing x, we see from Fig. 3, that this shadowed AKMS prediction increases more 
slowly than the unshadowed one, but that it keeps the characteristic "BFKL Q 2 curvature" . 

To conclude, we have performed a detailed analysis of the Q 2 dependence of the structure 
function F 2 (x,Q 2 ) in the small x region which is being probed at HERA. We have found 
that the theoretically-motivated BFKL-based predictions do indeed lead, in the HERA small 
x regime, to a more pronounced curvature of F 2 (x,Q 2 ) than those based on next-to-leading 
order Altarelli-Parisi evolution. The difference is illustrated in Fig. 3 by the comparison of the 
AKMS curve with that for MRS(H). From the figure we see that data at the smallest possible 
x values will be the most revealing. The measurements shown are from the 1992 run, but data 
with much higher luminosity, and at smaller x, will become available in the near future. Clearly 
the experimental identification of the characteristic BFKL Q 2 behaviour will pose a difficult, 
though hopefully not an impossible, task. 



3 To be precise, we take F 2 (background) = F 2 (xq = 0.1,Q 2 )(x/xo) 08 a form which is motivated by 
"soft" Pomeron Regge behaviour |L6) . Other reasonable choices of the background do not change our conclusions. 



5 



Acknowledgements 

We thank Dick Roberts and James Stirling for valuable discussions. Two of us (JK, ADM) 
thank the European Community for Fellowships and three of us (AJA, KGB, PJS) thank the 
Polish KBN - British Council collaborative research programme for partial support. This work 
has also been supported in part by Polish KBN grant no. 2 0198 91 01 and by the UK Science 
and Engineering Research Council. 

References 

[1] HI collaboration: I. Abt et al., preprint DESY-93-117, August 1993 

[2] ZEUS collaboration: M. Derrick et al., Phys. Lett. B316 (1993) 412 

[3] E.A. Kuraev, L.N. Lipatov and V.S. Fadin, Sov. Phys. JETP 45 (1977) 199; Ya.Ya. 
Balitsky and L.N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822 

[4] A.D. Martin, W.J. Stirling and R.G. Roberts, Phys. Lett. B306 (1993) 145 

[5] A.D. Martin, W.J. Stirling and R.G. Roberts, preprint RAL-93-077 (1993) 

[6] M. Gliick, E. Reya and A. Vogt, Z. Phys. C53 (1992) 127 

[7] A.J. Askew, J. Kwiecihski, A.D. Martin and P.J. Sutton, Phys. Rev D47 (1993) 3775 

[8] S. Catani, M. Ciafaloni and F. Hautmann, Phys. Lett. B242 (1990) 91; Nucl. Phys. B366 
(1991) 135 

[9] T. Jaroszewicz, Phys. Lett. B116 (1982) 291 

[10] J. Kwiecihski, Z.Phys. C29 (1985) 561; J.C. Collins and J. Kwiecihski, Nucl. Phys. B316 
(1989) 307 

[11] J. Kwiecihski, A.D. Martin and P.J. Sutton, Phys. Rev. D44 (1991) 2640 

[12] A.J. Askew, J. Kwiecihski, A.D. Martin and P.J. Sutton, Mod. Phys. Lett. A (in press); 
see also Durham preprint DTP/93/28, September 1993 

[13] J.R. Forshaw, Phys. Lett. B314 (1993) 425 

[14] HI collaboration: I. Abt et al., preprint DESY-93-146, October 1993 

[15] CTEQ collaboration: J. Botts et al., preprint MSU-HEP 93/18 (1993) 

[16] A. Donnachie and PV. Landshoff, Phys. Lett. B296 (1992) 227; Manchester preprint 
M/C-TH 93/11 (1993) 

[17] L.V. Gribov, E.M. Levin and M.G. Ryskin, Phys. Rep. 100 (1983) 1 

[18] G. Marchesini, Proc. of Workshop "QCD at 200TeV", Erice, 1990, eds. L. Cifarelli and 
Yu.L. Dokshitzer, Plenum Press, 1992, p. 183, and references therein 



6 



Figure Captions 



Fig. 1: The measurements of F 2 (x,Q 2 ) at Q 2 = 15 and 30 GeV 2 by the H10 and ZEUS0 col- 
laborations shown by closed and open data points respectively, with the statistical and 
systematic errors added in quadrature; the HI and ZEUS data have a global normal- 
ization uncertainty of ±8% and ±7% respectively. The continuous, dotted and dashed 
curves respectively correspond to the values of F 2 obtained from MRS(H)[|], GRV|6|] and 
MRS(D'_)0] partons. The curves that are shown as a sequence of small squares (trian- 
gles) correspond to the unshadowed (strong or "hot-spot" shadowing) AKMS predictions 
obtained by computing F 2 = f ® F 2 ^ + F 2 (background) as in ref.[]l2f and as described in 
the text. 

Fig. 2: Diagrammatic display of the fc^-factorization formula (2), which is symbolically of the 
form F 2 = f £g> F 2 °\ where / denotes the gluon ladder and F 2 ^ the quark box (and 
crossed box) amplitude. 

Fig. 3: The Q 2 dependence of F 2 (x, Q 2 ) at small x (note the shifts of scale between the plots 
at the different x values, which have been introduced for clarity). The curves are as in 
Fig. 1. Also shown are the measurements of the 1992 HERA run obtained by the ZEUS 
collaboration pfl (open points) and, by the HI collaboration [H] using their "electron" 



analysis (closed points). Only statistical errors of the data are shown. The ZEUS points 
shown on the x=0. 00098 curves are measured at an average x=0. 00085. A challenge for 
future experiments is to distinguish between curves like AKMS and MRS(H), both of 
which give a satisfactory description of the existing data. 
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